Weak magnetism and non-Fermi liquids near heavy-fermion critical points 
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This paper is concerned with the weak-moment magnetism in heavy-fermion materials and its 
relation to the non-Fermi liquid physics observed near the transition to the Fermi liquid. We explore 
the hypothesis that the primary fluctuations responsible for the non-Fermi liquid physics are those 
associated with the destruction of the large Fermi surface of the Fermi liquid. Magnetism is suggested 
to be a low-energy instability of the resulting small Fermi surface state. A concrete realization of 
this picture is provided by a fractionalized Fermi liquid state which has a small Fermi surface of 
conduction electrons, but also has other exotic excitations with interactions described by a gauge 
theory in its deconfined phase. Of particular interest is a three-dimensional fractionalized Fermi 
liquid with a spinon Fermi surface and a U(l) gauge structure. A direct second-order transition 
from this state to the conventional Fermi liquid is possible and involves a jump in the electron 
Fermi surface volume. The critical point displays non-Fermi liquid behavior. A magnetic phase 
may develop from a spin density wave instability of the spinon Fermi surface. This exotic magnetic 
metal may have a weak ordered moment although the local moments do not participate in the 
Fermi surface. Experimental signatures of this phase and implications for heavy-fermion systems 
are discussed. 



I. INTRODUCTION 

The competition between the Kondo effect and inter- 
moment exchange determines the physics of a large class 
of materials which have localized magnetic moments cou- 
pled to a separate set of conduction electron When 
the Kondo effect dominates, the low-energy physics is 
well described by Fermi liquid theory (albeit with heavily 
renormalized quasiparticle masses) . In contrast when the 
inter-moment exchange dominates, ordered magnetism 
typically results. 

A remarkable experimental property of such magnetic 
states is that the magnetism is often very weak - the 
ordered moment per site is much smaller than the mi- 
croscopic local moment that actually occupies each site. 
The traditional explanation of this feature is that the 
magnetism arises out of imperfectly Kondo-screened lo- 
cal moments. In other words, the magnetism is to be 
viewed as a spin density wave that develops out of the 
parent heavy Fermi liquid state. We will henceforth de- 
note such a state as SDW. Clearly a SDW state may be 
a small moment magnet. 

A different kind of magnetic metallic state is also pos- 
sible in heavy-fermion materials where the moments or- 
der at relatively large energy scales, and simply do not 
participate in the Fermi surface of the metal. In such 
a situation, the saturation moment in the ordered state 
would naively be large, i.e., of order the atomic moment. 

Often the distinction between these two kinds of mag- 
netic states can be made sharply: the two Fermi surfaces 
in the two states may have different topologies (albeit, 



the same volume modulo the volume of the Brillouin zone 
of the ordered state), so that they cannot be smoothly 
connected to one another. 

In recent years, a number of experiments have un- 
earthed some fascinating phenomena near the zero tem- 
perature (T) quantum transition between the heavy- 
fermion liquid and the magnetic metal. In particu- 
lar, many experiments do not fit easily @, y, 0, 
into a description in terms of an effective Gaussian the- 
ory for the spin density wave fluctuations, renormalized 
self-consistently by quartic interactions 0, 0, H, EJ El • 
This theory makes certain predictions on deviations from 
Fermi liquid behavior as the heavy Fermi liquid state 
approaches magnetic ordering induced by the condensa- 
tion of the spin density wave mode; those predictions 
are, however, in disagreement with experimental findings. 
This conflict raises the possibility that the magnetic state 
being accessed is not in the first category discussed above: 
a SDW emerging from a heavy Fermi liquid. Rather, it 
may be the second kind of magnetic metal where the local 
moments do not participate at all in the Fermi surface. 
In other words, the experiments suggest that the Kondo 
effect (crucial in forming the Fermi liquid state) is itself 
suppressed on approaching the magnetic state. 

This proposal clearly raises several serious puzzles. 
How do we correctly describe the non-Fermi liquid 
physics near the transition? If this non-Fermi liquid be- 
havior is accompanied by the suppression of the Kondo 
effect, how do we reconcile it with the weak moments 
found in the magnetic state? The traditional explana- 
tion for the weak magnetism is apparently in conflict 
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with the picture that the Kondo effect and the resultant 
heavy Fermi liquid state are destroyed on approaching 
the magnetic state. In other words, the naive expecta- 
tion of a large saturation moment in a magnetic metal 
where the local moments do not participate in the Fermi 
surface must be revisited. 

The weakness of the ordered moment in the magnetic 
state may be reconciled with the apparent suppression of 
the Kondo effect if we assume that there are strong quan- 
tum fluctuations of the spins that reduce their moment. 
Such strong quantum effects may appear to be unusual 
in three-dimensional systems, but may be facilitated by 
the coupling to the conduction electrons (even if there 
is no actual Kondo screening). In this paper we study 
specific states where such quantum fluctuations have sig- 
nificantly reduced the ordered moment (or even caused it 
to vanish), and the evolution of such states to the heavy 
Fermi liquid. 

We begin with several general pertinent observations. 
First, consider the heavy Fermi liquid state. This Fermi 
liquid behavior is accompanied by a Fermi surface which, 
remarkably, satisfies Luttinger's theorem only if the lo- 
cal moments are included as part of the electron count. 
(Such a Fermi surface is often referred to as the "large 
Fermi surface" , and we will henceforth refer to such 
a phase as FL). The absorption of the local moments 
into the Fermi volume is the lattice manifestation of the 
Kondo screening of the moments. We take as our starting 
point the assumption that the Kondo effect becomes sup- 
pressed on approaching the magnetic state. What then 
happens to the large Fermi surface? 

In thinking about the resulting state theoretically, it 
is important to realize that once magnetic order sets in, 
there is no sharp distinction between a large Fermi vol- 
ume which includes the local moments, and a Fermi vol- 
ume that excludes the local moments - the latter is often 
loosely referred to as "small" . This is because the Fermi 
volumes can only be defined modulo the volume of the 
Brillouin zone, and the onset of magnetic order at least 
doubles the unit cell and hence at least halves the Bril- 
louin zone volume. (There can, however, be a distinction 
between the Fermi surfaces topologies in the two situa- 
tions.) 

In this paper we will take the point of view that the 
primary transition involves the destruction of the large 
Fermi surface, and that the resulting small Fermi surface 
state has a distinct physical meaning even in the absence 
of magnetic order. The magnetic order will be viewed 
as a low-energy instability of the resulting state in which 
the local moments are not to be included in the Fermi 
volume. 

Evidence in support of this point of view exists. In 
the experiments the non-Fermi liquid behavior extends 
to temperatures well above the Neel ordering tempera- 
ture even far away from the critical point. This suggests 
that the fluctuations responsible for the non-Fermi liq- 
uid behavior have very little to do with the fluctuations 
of the magnetic order parameter. Some further support 



is provided by the results of inelastic neutron scatter- 
ing experiments that apparently see critical behavior at 
a range of wave-vectors including (but not restricted to) 
the one associated with magnetic ordering in the mag- 
netic metal 5]. Finally, there even exist materials in 
which the non-Fermi liquid features persist into the mag- 
netically ordered side - this is difficult to understand if 
the non-Fermi liquid physics is attributed to critical fluc- 
tuations of the magnetic order parameter. 

Conceptually, as we asserted above, it pays to allow for 
the possibility of a non-magnetic state in which the sup- 
pression of the Kondo effect removes the local moments 
from the Fermi volume, resulting in a "small Fermi sur- 
face" , even though such a state may not actually be a 
ground state in the system of interest. In our previous 
work we argued that such states do exist as ground 
states of Kondo lattice models in regular d-dimensional 
lattices, and that the violation of Luttinger's theorem in 
such a state was intimately linked to the presence of neu- 
tral S = 1/2 and 5 = excitations induced by topologi- 
cal order (see also Appendix |A"|) : we dubbed such ground 
states FL*. 

Clearly, it is worthwhile to explore metallic magnetic 
states that develop out of such FL* states (just as the 
usual SDW state develops out of the Fermi liquid) . Such 
states, which we will denote SDW*, represent a third 
class of metallic magnetic states distinct from both the 
conventional SDW and the conventional local-moment 
metal described above. As we will see, in such mag- 
netic states the local moments do not participate in the 
Fermi surface. Nevertheless they may have a weak or- 
dered moment. Thus these states offer an opportunity 
for resolution of the puzzles mentioned above. The prop- 
erties and the evolution of such states, and their parent 
FL* states, to the Fermi liquid will be the subject of this 
paper. The SDW* states inherit neutral spin S = 1/2 
spinon excitations and S — "gauge" excitations from 
the FL* states, which will be described more precisely 
below; these excitations coexist with the magnetism and 
the metallic behavior. The experimental distinction be- 
tween the SDW and SDW* states is however subtle, and 
will also be described in this paper. (The FL and FL* 
states can be easily distinguished by the volumes of the 
Fermi surfaces, but this distinction does not extend to 
the SDW and SDW* states.) 

We emphasize that a wide variety of heavy-fermion ma- 
terials display non-Fermi liquid physics in the vicinity of 
the onset of magnetism that is, to a considerable extent, 
universal. However, the detailed behavior at very low 
temperature appears to vary across different systems. In 
particular, in some materials a direct transition to the 
magnetic state at very low temperature does not occur 
(due for instance to intervention of a superconducting 
state). In other materials, such a direct transition does 
seem to occur at currently accessible temperatures. In 
view of this, we will not attempt to predict the detailed 
phase diagram at ultra-low temperatures. We focus in- 
stead on understanding the universal non-Fermi liquid 



3 



physics not too close to the transition and its relation to 
the magnetic state. 



A. Summary of results 

Our analysis is based upon non-magnetic translation- 
invariant states that have a small Fermi surface (FL*), 
and the related transitions to the heavy Fermi liquid 
(FL). As we showed previously ^lj, the FL* state has 
a Fermi surface of long-lived electron-like quasiparticles 
whose volume does not count the local moments. The 
local moments are instead in a state adiabatically con- 
nected to a spin-liquid state with emergent gauge exci- 
tations. Such spin liquids can be classified by the gauge 
group determining the quantum numbers carried by the 
neutral S — 1/2 spinon excitations and the gauge exci- 
tations, and previous work 0, has shown that the 
most prominent examples are Z 2 and U(l) spin liquids. 
The Z 2 spin liquids are stable in all spatial dimensions 
d > 2, while the U(l) spin liquids exist only in d > 3 (the 
latter correspond to the existence of a Coulomb phase in 
a compact U(l) gauge theory in d > 3, as discussed in 
Ref . fl4f) . Correspondingly, we also have the metallic Z 2 
FL* and U(l) FL* states. Our previous work [llj consid- 
ered primarily the Z 2 FL* state, whereas here we focus 
on the U(l) FL* state. 

As we have already discussed, these non-magnetic 
states may lead to magnetic order at low energies, or in 
proximate states in a generalized phase diagram. In this 
manner the FL state leads to the SDW state, while the 
FL* states lead to the Z 2 SDW* and the U(l) SDW* 
states. The relation between the metallic SDW and 
SDW* states has a parallel to that between the insulating 
Neel state and the AF* state of Refe.IT3.ll5l 

We will also discuss the evolution from the U(l) SDW* 
state to the conventional Fermi liquid. As explained ear- 
lier, the underlying transition is that between FL and 
FL* states which controls the nature of the Fermi sur- 
face. In Ref. IllL we argued that the spinon pairing in 
the Z 2 FL* state implied that there must be a super- 
conducting state in between the FL and Z 2 FL* states. 
There is no such pairing in the U(l) FL* state, and hence 
there is the possibility of a direct transition between the 
FL and U(l) FL* states: this transition and the nature 
of the states flanking it are the foci of our paper. Note 
that the volume of the Fermi surface jumps at this tran- 
sition. Nevertheless the transition may be second order. 
This is made possible by the vanishing of the quasipar- 
ticle residue on an entire portion of the Fermi surface (a 
"hot" Fermi surface) on approaching the transition from 
the FL side. Non-Fermi liquid physics is clearly to be 
expected at such a second order Fermi-volume changing 
transition. We reiterate that the U(l) FL* state is only 
believed to exist in d > 2. 

We study the FL and U(l) FL* states by the "slave" 
boson method, introduced in the context of the single- 
moment Kondo problem |l6j |. In this method, the con- 
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FIG. 1: Crossover phase diagram for the vicinity of the d = 3 
quantum transition involving breakdown of Kondo screening. 
Jk is the Kondo exchange in the Hamiltonian introduced in 
Section ITTTI The only true phase transition above is that at 
the T — quantum critical point at Jk = Jkc between the 
FL and FL* phases. The "slave" boson b measures the mix- 
ing between the local moments and the conduction electrons 
and is also described in Section ITTTI The crossovers are sim- 
ilar to those of a dilute Bose gas as a function of chemical 
potential and temperature, as discussed in Refs. ITsL 12(3 — the 
horizontal axis is a measure of the boson chemical potential 
fib- The boson is coupled to a compact U(l) gauge field; at 
T — this gauge field is in the Higgs/confining phase in the 
FL state, and in the deconfining/Coulomb phase in the FL* 
state. There is no phase transition at T > between a phase 
with (6) 7^ and a phase with (6) = because such a transi- 
tion is absent in a theory with a compact U(l) gauge field in 
d = 3 pH (the mean-field theories of Sections. Illll and ll V Cl do 
show such transitions, but these will turn into crossovers upon 
including gauge fluctuations) . The compactness of the gauge 
field therefore plays a role in the crossovers in the "renormal- 
ized classical" regime above the FL state (this has not been 
worked out in any detail here). However, the compactness 
is not expected to be crucial in the quantum-critical regime. 
The crossover line displayed between the FL and quantum 
critical regimes can be associated with the "coherence" tem- 
perature of the heavy Fermi liquid. At low T, as discussed in 
the text, there are likely to be additional phases associated 
with magnetic order (the SDW and SDW* phases) , and these 
are not shown above but are shown in Fig. HI they also appear 
in the mean-field phase diagram in Fig. [I] 



densation of the slave boson marks the onset of Kondo co- 
herence that characterizes the FL phase. In contrast the 
slave boson is not condensed in the FL* phase. Fluctua- 
tions about this mean-field description lead to the critical 
theory of the transition involving a propagating boson 
coupled to a compact U(l) gauge field, in the presence of 
damping from fermionic excitations. 

We note that earlier studies [l7l[T8j | of single- impurity 
problems found a temperature-induced mean-field tran- 
sition between a state in which the slave boson is con- 
densed (and hence the local moment is Kondo screened) 
and a state in which the boson has no condensate: how- 
ever, it was correctly argued that this transition is an 
artifact of the mean-field theory, and no sharp transition 
exists in the single-moment Kondo problem at T > 0. 
If we now naively generalize this single-impurity model 
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FIG. 2: Expected phase diagram and crossovers for the evolu- 
tion from the U(l) SDW* phase to the conventional FL. Two 
different transitions are genencally possible at zero tempera- 
ture: Upon moving from the SDW* towards the Fermi liquid, 
the fractionalization is lost first followed by the disappear- 
ance of magnetic order. Nevertheless the higher temperature 
behavior in the region marked 'quantum critical' is non-fermi 
liquid like, and controlled by the Fermi volume changing tran- 
sition from FL to FL* . This may be loosely associated to the 
breakdown of Kondo screening. 



to the lattice, we will find that the T = ground state 
always has Kondo screening. It is only upon including 
frustrating inter-moment exchange interactions - equiva- 
lent to having "dispersing" spinons - that it is possible to 
breakdown Kondo screening and reach a state in which 
the slave boson is not condensed. This transition is not 
an artifact of mean-field theory, we show here that it re- 
mains sharply defined in d = 3. 

Our analysis of the above d = 3 U(l) gauge theory 
leads to the schematic crossover phase diagram as a func- 
tion of the Kondo exchange Jk and T shown in FigGJ 

The crossover phase diagram in Fig. ^is similar to that 
of a dilute Bose function of chemical potential 

and temperature 19, 20] . Here the bosons are coupled to 
a U(l) gauge field, and this is important for many of the 
critical properties to be described in the body of the pa- 
per. Notably, in Fig. ^ the density of bosons is not fixed, 
and varies as a function of T, Jk and other couplings in 
the Hamiltonian. Indeed, the contours of constant boson 
density have a complicated structure, which are similar 
to those in Ref. [23. This variation in the boson density 
is a crucial distinction from earlier analyses [22l |23| of 
boson models coupled to damped U(l) gauge fields: in 
these earlier works, the boson density was fixed at a T- 
independent value. As we will see, allowing the boson 
density to vary changes the critical properties, and has 
significant consequences for the structure of the crossover 
phase diagram and for the T dependence of observables. 

We will show that non-Fermi liquid physics obtains in 
the quantum critical region of this transition. Further- 
more, we argue that fluctuation effects may lead to a spin 



density wave developing out of the spinon Fermi surface 
of the U(l) FL* phase, thereby obtaining the U(l) SDW* 
phase. The expected phase diagram and crossovers for 
the evolution from the U(l) SDW* phase to the FL phase 
is shown in Fig. |2 We examine few different kinds of 
such U(l) SDW* phases depending on the details of the 
spinon Fermi surface. We also describe a number of spe- 
cific experimental signatures of the U(l) SDW* phase 
which may help to distinguish it from more conventional 
magnetic metals. 



B. Relation to earlier work 

We have already mentioned a number of precursors to 
our ideas in our discussion so far. Here, we complete 
this by noting some other related developments in the 
literature. 

Early on, Andrei and Coleman [24| and Kagan et al. 
[2l| discussed the possibility of the decoupling of local 
moments and conduction electrons in Kondo lattice mod- 
els. Andrei and Coleman had the local moments in a 
spin-liquid state which is unstable to U(l) gauge fluctu- 
ations, and did not notice violation of Luttinger's theo- 
rem. The possibility of small electronic Fermi surfaces 
was noted by Kagan et al., but no connection was made 
to the requirement this imposes on emergent gauge exci- 
tations pT| . 

More recently, Burdin et al. described many aspects 
of the physics we are interested in a dynamical mean- 
field theory of a random Kondo lattice |2f| . In this work, 
they obtained a state in which local moments formed a 
spin liquid and stayed essentially decoupled from the con- 
duction electrons. They emphasized that the transition 
between such a state (which is the analog of our FL* 
states) and a conventional heavy Fermi liquid (the FL 
state) should be understood as a Fermi volume changing 
transition. However questions of emergent gauge struc- 
ture were not addressed by them. 

Demler et al. [2?J discussed fractionalized phases of 
Kondo lattice models. However, they did not consider 
any states with long-lived electron-like quasiparticles, as 
are present in the FL* phase. 

Recently Essler and Tsvelik [2^] discussed the fate of 
one-dimensional Mott insulator under a particular long- 
range inter-chain hopping. At intermediate tempera- 
tures, they obtain a state with a small Fermi surface, in 
that the Fermi surface volume does not count the local 
moments [29| . However, their construction does not lead 
to a state with emergent gauge excitations in higher di- 
mensions, and as they conclude, their state is unstable to 
magnetic order at low temperatures. We believe this low- 
T state is an ordinary SDW state, and any realizations 
of small Fermi surfaces at intermediate temperatures are 
remnants of one-dimensional physics. In contrast, all our 
constructions are genuinely higher-dimensional, and only 
work for d > 2. 

The physics of the destruction of the large Fermi sur- 
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face by the vanishing of Kondo screening has been ad- 
dressed in interesting recent works [30, HH using an 
"extended dynamical mean-field theory". We have ar- 
gued in our discussion above that vanishing of Kondo 
screening is conceptually quite a different transition from 
the onset of magnetic order; consistent with this expec- 
tation, Sun and Kotliar [3ll ] found two distinct points 
associated with these transitions. It is our contention 
that the critical theory of the FL to U(l) FL* transition 
(discussed in the present paper) is the d — 3 realization of 
the large-dimensional critical point with vanishing Kondo 
screening found by Sun and Kotliar. 



C. Outline 

The rest of the paper is organized as follows: In Sec- 
tion [HJ we briefly review the properties of various frac- 
tionalized Fermi liquids (FL*). A specific U(l) FL* state 
where the spinons form a Fermi surface is considered. In 
Section IlIII we construct a mean-field description of this 
state and its transition to the heavy Fermi liquid. This 
transition involves a jump in the Fermi surface volume 
but is nevertheless shown to be second order within the 
mean-field theory. This is made possible by the vanishing 
of the quasiparticle residue Z on an entire Fermi surface 
(a "hot" Fermi surface) as one moves from the heavy 
Fermi liquid to the fractionalized Fermi liquid. Fluctua- 
tions about this mean-field description are then consid- 
ered. In Section ITVl we first consider fluctuation effects 
on the phases - in particular the FL* phase. We ar- 
gue that the specific heat coefficient 7 diverges logarith- 
mically once the leading-order fluctuations are included. 
Furthermore, fluctuations also make possible a spin den- 
sity wave instability of the spinon Fermi surface, leading 
to a U(l) SDW* state. To illustrate possible phases, we 
will discuss an improved mean-field theory which includes 
the SDW order parameter, and present phase diagrams 
showing the influence of temperature and magnetic field. 
We then examine fluctuation effects at the critical point 
of the transition between FL* and FL in Section [V] We 
argue that the logarithmic divergence of the specific heat 
coefficient persists in the quantum critical region, and 
also that non-Fermi liquid transport obtains there. In 
Section IVT1 we discuss the properties of the U(l) SDW* 
phase in greater detail with particular attention to its 
identification in experiments. A discussion of the impli- 
cations for various experiments in Sec. I VI il will conclude 
the paper. 



II. FRACTIONALIZED FERMI LIQUIDS 

The existence of non-magnetic translation invariant 
"small Fermi surface" states was shown in a recent article 
by us with a focus on two-dimensional Kondo lat- 
tices. Such states were obtained when the local-moment 
system settles into a fractionalized spin liquid (rather 



than a magnetically ordered state) due to inter-moment 
interactions. A weak Kondo coupling to conduction elec- 
trons does not disrupt the structure of the spin liquid 
but leaves a sharp Fermi surface of quasiparticles whose 
volume counts the conduction electron density alone (a 
small Fermi surface). Thus these states have fraction- 
alized excitations that coexist with conventional Fermi- 
liquid- like quasiparticle excitations. We dubbed these 
states FL* (to distinguish them from the conventional 
Fermi liquid FL). We also pointed out an intimate con- 
nection between the disappearance of the large Fermi sur- 
face and fractionalization, and this is discussed further in 
Appendix 1X1 

The FL* phase can be further classified by the nature 
of the spin liquid formed by the local moments. Recent 
years have seen considerable progress in our understand- 
ing of fractionalized spin liquids. An important feature of 
spin liquid states in d > 2 is that they possess emergent 
gauge structure. Put simply, this means that the dis- 
tinct excitations in such phases interact with each other 
through long ranged interactions which can be mathe- 
matically encapsulated as gauge interactions. In other 
words, the effective field theory of the state is a gauge 
theory in its deconfined phase. The two natural possi- 
bilities are that the emergent gauge group is either Zi 
or U(l). The former is allowed in any dimension d > 2 
while the latter is only allowed in d = 3 (or higher) . 

The Z2 states have been discussed at length in the lit- 
erature and in the present context in our earlier work 
[TT| . In contrast, the U(l) states have not been dis- 
cussed much, though their possible occurrence (in d = 3) 
and their universal properties have been appreciated by 
many workers in the field. We therefore provide a quick 
discussion: The distinct excitations in the d = 3 U(l) 
spin liquid phases are neutral spin- 1/2 spinons, a gap- 
less (emergent) gauge photon, and a gapped point defect 
(the "monopole"). The spinons are minimally coupled 
to the photon and hence interact through emergent long 
ranged interactions. For simple microscopic models that 
realize such phases, see Ref . fl4L l33l A crucial distinction 
between the Zi spin liquids is that the spinons in this 
phase are not generically paired, i.e., the spinon number 
is conserved |34j . 

Several classes of spin liquids are theoretically possible 
with the same gauge structure. These may be charac- 
terized by the statistics of the spinons, their band struc- 
ture, etc. For the rest of this paper, we will focus on a 
particular three-dimensional U(l) spin liquid state with 
fermionic spinons that form a Fermi surface. A specific 
toy model which displays this phase is presented in Ap- 
pendix [5] 

As with the Z2 spin liquids discussed in Ref. 0, the 
gauge structure in the U(l) spin liquid state is also stable 
to a weak Kondo coupling to conduction electrons [35| . 
The resulting U(l) FL* state consists of a spinon Fermi 
surface coexisting with a separate Fermi surface of con- 
duction electrons. There will also be gapless photon 
and gapped monopole excitations. The physical electron 



6 



Fermi surface (as measured by de-Haas van Alphen ex- 
periments for instance) will have a small volume that is 
determined by the conduction electrons alone. 

In our previous work, we pointed out that the transi- 
tion from a Z2 FL* phase to the heavy FL will generically 
be preempted by superconductivity. This is due to the 
pairing of spinons in the Z2 phase. In contrast, we ex- 
pect that due to conservation of spinon number a direct 
transition between the U(l) FL* and heavy FL phases 
should be possible. 



III. MEAN-FIELD THEORY 

A simple mean-field theory allows a description both 
of a U(l) FL* phase and its transition to the heavy FL. 
Consider a three-dimensional Kondo-Heisenberg model, 
for concreteness on a cubic lattice: 



k r 

+ Jh y S r ■ S r ' . 

{rr') 



olol' War 



(1) 



Here c ka represent the conduction electrons and S r the 
spin-1/2 local moments on the sites of a cubic lattice, 
summation over repeated spin indices a is implicit. We 
use a fermionic "slave-particle" representation of the lo- 
cal moments: 



1 



Sr — nfra^aa'fr 



(2) 



where f ra describes a spinful fermion destruction opera- 
tor at site r. 

Proceeding as usual, we consider a decoupling of both 
the Kondo and Heisenberg exchange using two auxiliary 
fields in the particle-hole channel. Treating the fluctua- 
tions of these auxiliary fields by a saddle point approxi- 
mation (formally justified for a large- N SU(iV) general- 
ization), we obtain the mean-field Hamiltonian 



, tkC' ka C ka - XO (frafr'a + h.C. 
k (rr') 



M/ tijra ~ b ^ (c\j ka + h.C.) (3) 



where we assumed xo and b to be real, and have dropped 
additional constants to H. The mean-field parameters 
&o,X0 7M/ are determined by the conditions 



1 — ifrafra) , 
2&0 = Jl<(ci a fra}, 
2X0 = Jniflafr'a)- 



(4) 
(5) 
(6) 



In the last equation r, r' are nearest neighbors. 

There are two qualitatively different zero-temperature 
phases. First, there is the usual Fermi liquid (FL) phase 



when 6o,XO:A*/ are all non-zero. (Note that 60 ^ im- 
plies that Xo 7^ 0). This phase is readily seen to have a 
large Fermi surface as expected. Second, there is a phase 
(FL*) where the Kondo hybridization &o = but xo 7^ 0. 
(In this phase /i/ = 0.) This mean-field state repre- 
sents a situation where the conduction electrons are de- 
coupled from the local moments and form a small Fermi 
surface. The local-moment system is described as a spin 
fluid with a Fermi surface of neutral spinons. We expect 
that xo ~ Jh- 

The transition between these two different states can 
also be examined within the mean-field theory. Interest- 
ingly, the transition is second order (despite the jump in 
Fermi volume) and is described by 60 - * on approach- 
ing it from the Fermi liquid side. How can a second order 
transition be associated with a jump in the volume of 
the electron Fermi surface? This can be understood by 
examining the Fermi surfaces closely in this mean-field 
theory. 

The mean-field Hamiltonian is diagonalized by the 
transformation 



Cka = U k J ka+ + V k J ka . 
fka — V k "f ka+ — U k J ka . 



(7) 



Here "f ka ± are new fermionic operators in terms of which 
the Hamiltonian takes the form 

Hmi = Y E k+ll a + lk a + + £'fe-7| a _7fea- , (8) 



with 



E k± = 



e k + e k f ± f e k -e kf 



+ bl (9) 



Here e kf = fif - Xo J2 a =i,2,3 cos ( k a)- The u k ,v k intro- 
duced above are determined by 



u k 



bov k 



E k+ — e k 



(10) 



Consider first the FL* phase where 60 = = but 
Xo 7^ 0. The electron Fermi surface is determined by 
the conduction electron dispersion e k and is small. The 
spinon Fermi surface encloses one spinon per site and has 
volume half that of the Brillouin zone. For concreteness, 
we will consider the situation where the electron Fermi 
surface does not intersect the spinon Fermi surface. We 
will also assume that the conduction electron filling is 
less than half. 

Now consider the FL phase near the transition (small 
60). In this case, there are two bands corresponding to 
E k ± : one derives from the c-electrons (with weak / char- 
acter) while the other derives from the /-particles (with 
weak c-character) . We will call the former the c-band and 
the latter the /-band. For small bo, both bands intersect 
the Fermi energy so that the Fermi surface consists of 
two sheets (see Fig. EJ>- The total volume is large, i.e, 
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This then gives 



Spinon Fermi 
surface 



Small electron 
fermi surface 



Fermi Liquid 



Fractionalized Fermi Liquid 



FIG. 3: Fermi surface evolution from FL to FL*: close to the 
transition, the FL phase features two Fermi surface sheets 
(the cold c and the hot / sheet, see text). Upon approach- 
ing the transition, the quasiparticle residue Z on the hot / 
sheet vanishes. On the FL* side, the / sheet becomes the 
spinon Fermi surface, whereas the c sheet is simply the small 
conduction electron Fermi surface. 



includes both local moments and conduction electrons. 
Upon moving toward the transition to FL* (60 decreas- 
ing to zero) , the c- Fermi surface expands in size to match 
onto the small Fermi surface of FL* . On the other hand, 
the /-Fermi surface shrinks to match onto the spinon 
Fermi surface of FL*. 

Upon increasing 60 m the FL state and depending on 
the band structure, another transition is possible, where 
the c band becomes completely empty. Then, the Fermi 
surface topology changes from two sheets to a single sheet 
- such a transition between two conventional Fermi liq- 
uids is known as Lifshitz transition and will not be further 
considered here. 

The quasiparticle weight Z close to the FL-FL* transi- 
tion is readily calculated in the present mean-field theory. 
For the electron Green's function we find 



-E. 



k+ 



— Ek 



(11) 



Therefore at the Fermi surface of the c-band (which has 
dispersion Ek+, the quasiparticle residue Z = u\. At this 
Fermi surface, Ek+ ~ ek » so that 



Ek+ ~ €k + 



£fe - ekf 



Uk 



H 



b 



(12) 



Using Eqs. 1)100. we then find Z as 1 on the c- Fermi 
surface. 

At the Fermi surface of the /-band on the other hand, 
Z = v\. Also near this Fermi surface, \ek — €kf\ ~ t 
where t is the conduction electron bandwidth. We have 
assumed as is reasonable that t ^S> Jh- Thus for the 
/-Fermi surface, 



Ek+ ~ £fc 



f 

T 



Vk- 



(13) 



vl 



(14) 



Thus the quasiparticle residue stays non-zero on the 
Fermi surface while it decreases continuously to zero on 
the /-Fermi surface on moving from FL to FL*. (The 
/-Fermi surface is "hot" while the c-Fermi surface is 
"cold".) 

Clearly the critical point is not a Fermi liquid. Z van- 
ishes throughout the hot Fermi surface at the transition, 
and non- Fermi liquid behavior results. It is interesting 
to contrast this result with the spin-fluctuation model 
(Hertz- Moriya-Millis criticality) where the non- Fermi liq- 
uid behavior is only associated with some "hot" lines in 
the Fermi surface, and consequently plays a subdominant 
role. 

Despite the vanishing quasiparticle weight Z, the effec- 
tive mass m* of the large Fermi surface state does not di- 
verge at the transition in this mean-field calculation, be- 
cause the electron self-energy is momentum-dependent. 
Physically, the quasiparticle at the hot Fermi surface is 
essentially made up of the /-particle for small b; even 
when b goes to zero the /-particle (the spinon) continues 
to disperse due to the non-vanishing xo term. Indeed the 
low-temperature specific heat C ~ 7T with 7 non-zero 
in both phases. As we argue below, this is an artifact 
of the mean-field approximation and will be modified by 
fluctuations. 

The detailed shape of the spinon Fermi surface in the 
FL* phase (or the hot Fermi surface which derives from 
it in the FL phase) depends on the details of the lat- 
tice and the form of the local moment interactions. For 
the particular model discussed above, the spinon Fermi 
surface is perfectly nested. In more general situations, a 
non-nested spinon Fermi surface will obtain. In all cases, 
however, the volume of the spinon Fermi surface will cor- 
respond to one spinon per site. 



IV. FLUCTUATIONS: MAGNETISM AND 
SINGULAR SPECIFIC HEAT 

Fluctuation effects modify the picture obtained in the 
mean-field theory in several important ways. We first 
discuss fluctuation effects in the two phases. The heavy 
Fermi liquid phase is of course stable to fluctuations 
- their main effect being to endow the /-particle with 
a ph ysic al electric charge thereby making it an elec- 
tron [3a, |33 . Fluctuation effects are more interesting in 
the FL* state, and are described by a U(l) gauge theory 
minimally coupled to the spinon Fermi surface (which 
continues to be essentially decoupled from the conduc- 
tion electron small Fermi surface). This may be made 
explicit by parameterizing the fluctuations in the action 
in the FL* phase as follows: 



Xrr> (r) 



e la " AT) Xor 



(15) 
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The action then becomes 



k^ kj I k 



s = 
s c = 

Sf = 



S c f - 



Sb 



S c + Sf + Sf c + Si, , 

/ dT ^2 d k(d T - £fc)Cfe , 
^ k 

/ dT^/ r (5 r - ia )f r 
J r 

-E^ (e w ""7r/,.' + h.c.) 

(rr 1 ) 

~ j dT E ( b rC r fr + h.C.) , 



(16) 



T\w n \/k + X fk 2 



(17) 



d -E 



4|6 r 



Jk 



As usual, the field clq is introduced to impose the con- 
straint that there is one spinon per site and may be inter- 
preted as the time component of the gauge field. By as- 
sumption b r is not condensed. It is useful to start by com- 
pletely ignoring all coupling between c and / fermions. 
The action for the / particles describes a Fermi surface 
of spinons coupled to a compact U(l) gauge field. 

An important simplification for the three-dimensional 
systems of interest (as compared to d = 2) is that the 
U(l) gauge theory admits a deconfined phase where the 
spinons potentially survive as good excitations of the 
phase. In what follows we will assume that the system 
is in such a deconfined phase. (This is formally justified 
in the same large- TV limit as the one for the mean field 
approximation.) This deconfined phase has a Fermi sur- 
face of spinons coupled minimally to a gapless "photon" 
(U(l) gauge field). (Due to the compactness of the un- 
derlying gauge theory, there is also a gapped monopole 
excitation.) Thus two static spinons interact with each 
other through an emergent long range 1/r Coulomb in- 
teraction. Putting back a small coupling between the 
c and / particles will not change the deconfined nature 
of this phase. (In particular the monopole gap will be 
preserved.) This is the advocated U(l) FL* phase. 



A. Specific heat 

The coupling of the massless gauge photon to the 
spinon Fermi surface leads to several interesting modi- 
fications of the mean field results. First, consider the 
effect of the spatial components of the gauge field. It is 
useful to work in the gauge V ■ a = so that the vector 
potential is purely transverse. Unless otherwise stated, 
we assume a generic spinon Fermi surface (without flat 
portions) henceforth. Integrating out the spinons and ex- 
panding the resulting action to quadratic order gives the 
following well-known form for the propagator for these 
transverse gauge fluctuations: 

D i:j (k,iuj n ) = (ai(k,ico n )aj(-k,-iu} n )} 



Here T, \f are positive constants that are determined by 
the details of the spinon dispersion, and w„ is an imag- 
inary Matsubara frequency. Note that the gauge fluctu- 
ations are overdamped in the small q limit. As was first 
shown in a different context by Holstein et al. (and re- 
viewed in Appendix[DJ) , this form of the gauge field action 
leads to a T In 1 /T singularity in the low-temperature 
specific heat. Thus the specific heat coefficient 7 = C/T 
diverges logarithmically at low temperature in the U(l) 
FL* phase. 

We also briefly mention the effect of the longitudinal 
(time-component) of the gauge field. This couples to the 
local / fermion density, and so its influence is very much 
like a repulsive density-density interaction. The longitu- 
dinal gauge field propagator has a structure very similar 
to that of a standard RPA density fluctuation propaga- 
tor, and so does not lead to any non-Fermi liquid behav- 



B. Magnetic instability 

The repulsive interaction mediated by the longitudinal 
part of the gauge interaction can lead to various instabil- 
ities of the spinon Fermi surface. In particular, it is inter- 
esting to consider an SDW instability of the spinon Fermi 
surface. The resulting state will have magnetic long range 
order that could potentially have a weak moment as it is 
an SDW state that is formed out of the spinon Fermi sur- 
face. However, in contrast to the traditional view of the 
weak magnetism, here the SDW instability is not that of 
the large Fermi surface heavy Fermi liquid. Despite the 
occurrence of magnetic long range order, this magnetic 
state is far from conventional. Because the SDW order 
parameter is gauge neutral, the presence or absence of a 
SDW condensate has little substantive effect on the struc- 
ture of the gauge fluctuations. Indeed, the latter remain 
as in the U(l) FL* state, even after the magnetic order 
has appeared in the descendant U(l) SDW* state. The 
spinons continue to be deconfined and are coupled to a 
gapless U(l) gauge field. Further, the monopole survives 
as a gapped excitation - this yields a sharp distinction 
with more conventional magnetic phases. These gauge 
excitations coexist with the gapless magnons associated 
with broken spin rotation invariance and with a Fermi 
surface of the conduction electrons. However, due to the 
broken translational symmetry in this state, there is no 
sharp distinction between small and large Fermi surfaces. 
So to reiterate, the exotic magnetic metal, dubbed U(l) 
SDW* , emerges as a low-energy instability of the spinon 
Fermi surface of the parent U(l) FL* state. 

Different possibilities emerge for the formation of the 
spin density wave out of the parent U(l) FL* phase, de- 
pending on the details of the spinon Fermi surface and 
the strength of the interactions driving the SDW insta- 
bility. We enumerate some of them below: 
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• (A) Perfectly nested spinon Fermi surface: 

In this case, arbitrarily weak interactions will drive 
an SDW instability. In the resulting state, the 
spinons are gapped. So upon integrating out the 
spinons, the effective action for the gauge field can 
be expanded safely in spatial and temporal gradi- 
ents, with no long-range couplings. Gauge invari- 
ance now demands that these terms in the gauge 
field action have the standard Maxwell form. Con- 
sequently, the photon becomes a sharp propagating 
mode at low energies (below the spinon gap) with 
linear dispersion. Despite clearly being a distinct 
phase from conventional spin density wave metals, 
the experimental distinction is subtle. 

• (B) Generic spinon Fermi surface, weak interaction: 

For a generic spinon Fermi surface, the leading spin 
density wave instability (which will require an in- 
teraction strength beyond some threshold value) 
will be at a wavevector that matches one of the 
"2fcp" wavevectors of the spinon Fermi surface. In 
the resulting state, a portion of the spinon Fermi 
surface (away from points connected by the order- 
ing wavevector) survives intact. The damping of 
the gapless U(l) gauge fluctuations due to coupling 
to gapless spinons is preserved. Consequently the 
low-temperature specific heat will continue to be- 
have as C(T) ~ Tln(l/T). Thus for this particu- 
lar U(l) SDW* state its non-Fermi liquid nature is 
readily manifested by specific heat measurements, 
providing a concrete example of a weak moment 
SDW metal with non-Fermi liquid thermodynam- 
ics at low temperature. 

• (C) Generic spinon Fermi surface, strong interac- 
tion: 

If the interactions are strong enough, even for a 
non-nested spinon Fermi surface, the spinons can 
develop a full gap with no portion of their Fermi 
surface remaining intact. The resulting phase is 
the same as that obtained in (A), and has a sharp 
propagating linear dispersing photon at low ener- 
gies. 

In Section Ivll we discuss experimental probes that can 
help distinguish these U(l) SDW* phase from the con- 
ventional spin density wave metals. 

C. Mean-field theory with magnetism 

In view of the possible occurrence of SDW phases we 
will now consider a modified mean-field theory which cap- 
tures the magnetic instability at the mean-field level, but 
does no longer correspond to a large- AT saddle point. We 
will discuss the fully self-consistent solution of the mean- 
field equations for arbitrary temperature and external 
magnetic field. 
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FIG. 4: Mean-field phase diagram of H m { 1181 on the cubic 
lattice, as function of Kondo coupling Jk and temperature 
T. Parameter values are electron hopping t = 1, Heisenberg 
interaction Jh = 0.1, decoupling parameter x = 0.2, and 
conduction band filling n c = 0.7. Thin (thick) lines are sec- 
ond (first) order transitions. The "decoupled" phase is an 
artifact of the mean-field theory, and the corresponding tran- 
sitions will become crossovers upon including fluctuations, as 
will the transition between the FL and U(l) FL* phases; the 
transitions surrounding the SDW and SDW* phases will of 
course survive beyond mean-field theory. 
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FIG. 5: Staggered magnetization determined from the mean- 
field solution H m f 11811 . Parameter are as in Fig. |1J the two 
curves correspond to two horizontal cuts of the phase diagram 
in Fig. E] At T = 0, the first-order character of the SDW- 
FL transition is clearly seen. Note that smaller values of the 
decoupling parameter x yield smaller values of the magneti- 
zation in the SDW and SDW* phases. 



The mean-field Hamiltonian, written down explicitly 
for SU(2) symmetry, takes the following form: 

H mt = J2£kc{ a C ka ~ (Xrr'fLfr'a + h.C.) 

k (rr 1 ) 

+ Vf,rfrafra ~ ^2 K (4 a f ra + h.C.) 

r r 
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2 ^^{H c S,r + H cxt ) ■ fl a O aa ' fra' 
r 

C ra (T act f C TC n' ~h ^const 

(18) 



where -ff GXt is the external field, and we have allowed for 
a spatial dependence of the mean- field parameters /U/,r, 
Xrr'i b r , ffeff,r- They have to be determined from the 
following equations: 



1 ifrafra) j 

26 r = J K{cl a f ra) , 
2Xrr' = (1 - x)J H {fl a fr'a) 



(19) 

(20) 
(21) 



i?cff,r = ^H^M r , , M r = -{fl a Saaffraf) , (22) 



where the last sum runs over the nearest neighbors r' of 
site r. We have introduced a parameter a; which allows 
to control the balance between ordered local-moment 
magnetism and spin-liquid behavior of the / electrons. 
A value x = 1/2 would correspond to an unrestricted 
Hartree-Fock treatment of the original Heisenberg inter- 
action; we will employ values x < 1/2 in order to model 
a weak magnetic instability of the spinon Fermi surface 
state. The constant piece of the Hamiltonian reads 



E c 



E/^+E^ 



(23) 



I (1 - x)J H 

For simplicity, we consider a simple cubic lattice, and 
assume a tight-binding dispersion for the conduction elec- 
trons, efe = — 2i 2Zo,=i 2 3 cos (^a) — Mcj where /j, c controls 
the conduction band filling. The mean-field equations 
can be self-consistently solved using a large unit cell, al- 
lowing for spatially inhomogeneous phases |39| . In this 
section we restrict our attention to mean-field solutions 
where the Xrr' = Xo fields are real (time-reversal invari- 
ant) and obey the full lattice symmetries, and b r = bo is 
site-independent. We employ a 2 x 1 unit cell, then anti- 
ferromagnetism is characterized by M r -x = M s exp(iQ-r) 
where Q = (7r,7r,7r) is the antiferromagnetic wavevector, 
and x is the magnetization axis (which is arbitrary in 
zero external field). 

In Fig. 01 we show a phase diagram obtained from self- 
consistently solving (|18fl together with the above mean- 
field equations at zero external magnetic field. A U(l) 
FL* phase with bo — and xo 7^ is realized at interme- 
diate temperatures. As expected, it is unstable to mag- 
netic order at low T, resulting in a U(l) SDW* ground 
state for small Jjc - this phase has in addition M s ^ 0. 
For the present parameter values, the spinon Fermi sur- 
face is gapped out in the SDW* phase. Increasing Jk 
drives the system into the FL phase with bo ^ 0, xo ^ 0, 
and M s = 0; at low temperatures a conventional SDW 



phase intervenes where all bo, xo, M s are non-zero. Note 
that the transition between FL and SDW is weakly first 
order at low temperatures. At high temperature, the 
mean-field theory only has a "decoupled" solution with 
= Xo — M s = - this decoupling is a well-known 
mean-field artifact and reflects the presence of incoher- 
ent excitations. 

In the FL phase, the above mentioned Lifshitz tran- 
sition occurs at Jk ~ 1.7 in the low-temperature limit, 
i.e., for Jk > 1.7 only a single Fermi surface sheet re- 
mains. Note that this transition does not lead to strong 
singularities in the mean-field parameters. 

The staggered magnetization of the SDW and SDW* 
states as determined from the mean-field solution are 
shown in Fig. [5] we can expect that fluctuation correc- 
tions will significantly reduce these mean-field values. We 
have also studied different values of the decoupling pa- 
rameter x: in particular smaller values of x lead to a 
suppression of ordered magnetism in favor of the non- 
magnetic FL* state, i.e., the SDW instability of FL* is 
shifted to lower temperatures (and becomes completely 
suppressed at small x); similarly, the ordered moment in 
the SDW phases is decreased with decreasing x. 

Interesting physics obtains when an external magnetic 
field is turned on, and the corresponding mean-field phase 
diagram is discussed in Appendix IU1 



V. FLUCTUATIONS NEAR THE FERMI 
VOLUME CHANGING TRANSITION 

We now turn to the effects of fluctuations beyond the 
mean-field theory at the phase transition between the FL 
and U(l) FL* phases. In mean-field theory, this transi- 
tion occurs through the condensation of the slave boson 
field b. Such a condensation survives as a sharp transition 
beyond mean-field only when T = 0. 

We begin by observing that in the mean-field theory all 
the important changes near the transition occur at the 
hot Fermi surface. The cold Fermi surface (essentially 
made up of c-particles) plays a spectator role. We there- 
fore integrate out the c-fields completely from the action 
in Eq. i|16[l to obtain an effective action involving the 
b, f and gauge fields alone. We also partially integrate 
out / excitations well away from the hot Fermi surface: 
this changes the b effective action from the simple local 
term in lltj|) . and endows it with frequency and momen- 
tum dependence. In this manner we obtain the following 
effective action at long distance and time scales: 



S — Sb + Sf , 
.S',, = / drd : V 



(24) 



T ( B (Vr ~ idf , 

b I o T - fib - ia I o 



2m& 



';\b\ 4 



(25) 



and Sf has the same form as in (|16H . Notice that the 
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b field has become a propagating boson, with the same 
terms in the action as a microscopic canonical boson: 
here these terms arise from a (6,/) fermion polarization 
loop integrated well away from the / Fermi surface. The 
parameters fj,t,,mb may be interpreted as the chemical 
potential and mass of the bosons respectively. The (b,f) 
fermion loop will also lead to higher time and spatial gra- 
dient terms as well as a density-density coupling between 
b and / in l|24(l . but all these are formally irrelevant near 
the quantum critical point of interest. 

A key feature of Q24)). induced by taking the spatial 
and temporal continuum limit, is that we have lost in- 
formation on the compactness of the U(l) gauge field 
a, i.e., the continuum action is now no longer periodic 
under a rr i — > a rr > + 2ir, as was the case for the lattice 
action ifTBl) . The U(l) gauge field is now effectively non- 
compact, and consequently monopole excitations have 
been suppressed. The monopole gap is finite in the U(l) 
FL* phase (which is the analo g of the "Coulomb" phase 
of the compact gauge theory) [4(|. In the FL phase, the 
monopoles do not exist - they are confined to each other. 
This occurs due to the condensation of the boson field. 
However, the monopole gap is not expected to close at 
the transition [42], and so neglecting the compactness 
of the gauge field is permissible. Indeed, the continuum 
action (|24[) provides a satisfactory description of the crit- 
ical properties of the FL to U(l) FL* transition at T = 0. 
However, as we noted in the caption of Fig^ the com- 
pactness of the gauge field is crucial in understanding the 
absence of a T > phase transition above the FL phase 

m 

The action in Eq. I)24|) above is similar to that pop- 
ular in gauge theory descriptions [2^, H3 of the normal 
state of optimally doped cuprates but with some crucial 
differences. Here the chemical potential of the bosons is 
fixed while in Refs. l22T.l23l the boson density was fixed; as 
we will see below, this significantly modifies the physical 
implications of the critical theory, and the nature of the 
non- Fermi liquid critical singularities as T > 0. Further- 
more, we are interested specifically in d = 3, as opposed 
to the d = 2 case considered in Refs. I22II23I 

The phase diagram of the action (|24l) was sketched in 
FigD The horizontal axis, represented in Fig ^ by Jk, 
is now accessed by varying Without any additional 
(formally irrelevant) second-order time derivative terms 
for b in the action, the quantum critical point between 
the FL and U(l) FL* phases occurs precisely at /15 = 0, 
T = 0. We will now discuss the physical properties in 
the vicinity of this critical point first at T = 0, and then 
at T > 0, followed by an analysis of transport properties 
using the quantum Boltzmann equation in Section IV CI 
The final subsection IV Dl will comment on the effect of 
the SDW or SDW* phases that may appear at very low 
temperatures (these are not shown in Fig.^ but sketched 
in Fig.©. 



A. Zero temperature 

In a mean-field analysis of l|24|l. we see that the FL* 
phase (the "Coulomb" phase of the gauge theory) obtains 
for W, < with (b) = 0, while the FL phase (the "Higgs" 
phase of the gauge theory) obtains for fib > 0. 

Consider fluctuations for fib < in the FL* phase. 
Here, there are no bosons in the ground state, and all self- 
energy corrections associated with the quartic coupling u 
vanish 41] . The gauge field propagator is given by l|17fl . 
and this does contribute a non-zero boson self energy. At 
small momenta p and imaginary frequencies e, the boson 
self-energy has the structure (determined from a single 
gauge-boson exchange process, as in Refs. 12^, 113) 

E 6 (fc, ie) ~ fc 2 (l + Cl \e\ ln(l/|e]) + ...), (26) 

where c\ is some constant. Apart from terms which 
renormalize the boson mass mb, these self-energy correc- 
tions are less relevant than the bare terms in the action, 
and so can be safely neglected near the critical point. No- 
tice also that £&(0,0) = 0, and so the quantum critical 
point remains at fib = 0. 

The critical exponents can now be determined as in 

Refs. nam 

and are simply those of the mean-field theory 

of gQ: 

v = l/2 ; z = 2 ; 77 = 0. (27) 

As in Il26[l we can also determine the fate the boson quasi- 
particle pole as influenced by the gauge fluctuations; we 
obtain 

ImS 6 (k, e = ~ sgn(e)e 2 ln(l/|e|). (28) 

The boson lifetime is clearly longer than its energy, and 
this pole remains well defined. Finally, we recall our 
statement in Section IIV Al that the gauge fluctuations 
lead to a Tln(l/T) specific heat in the FL* phase, with 
a diverging 7 co-efficient. This behavior remains all the 
way up to, and including, the critical point. Parentheti- 
cally, we note that the same calculation in d = 2 dimen- 
sions will yield C oc T 2 / 3 . 

We turn next to fib > 0, in the FL phase. Here the 
bosons are condensed, and Q26[l or explicit calculations 
show that 

(b) = b ~ {fib) 1 ' 2 ~ (Jk - Jkc) 1 ' 2 , (29) 

where Jkc is the position of the critical point in Fig ^ 
The transverse gauge field propagator may be obtained 
as in Section llV Al bv integrating out both the bosons and 
fermions and expanding the resulting action to quadratic 
order; the boson condensate leads to a "Meissner" term 
in the gauge propagator so that l|17fl is replaced by 

D i3 (k,iLu n ) = (ai{k,iuj n )aj(—k,—u n )) 
Sjj — kjkj/k 2 
T\Lo n \/k + X fk2+p s - 1 ' 
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Here p s is the boson "superfiuid density", and we have 
p s ~ 6g. The presence of such a Meissner term cuts off 
the singular gauge fluctuations. The divergence of the 
specific heat coefficient j(T) as a function of temperature 
at the critical point implies that it diverges at T = on 
approaching the transition from the FL side. As shown 
in Appendix IdI this is indeed the case, and we find that 7 
diverges as 7 ~ ln(l/6o)- In experiments, such a diverg- 
ing 7 is sometimes interpreted as a diverging effective 
mass. Importantly, the divergence of 7 is unrelated to 
the singularity in the quasiparticle residue on the "hot" 
Fermi sheet, Z, which obeys Z ~ 6q as shown in 1)140. 
and so vanishes linearly as a function of Jk — Jkc- 



B. Non-zero temperatures 

A crucial change at T > is that it is now no-longer 
true that £(,(0,0) = in a region with (b) = 0. Instead, 
as in earlier studies of the dilute Bose gas 0, |2j| , we 
have 



£ fc (0,0) = 2u 



d d k 



1 



(27r) d exp [jfe 2 /(2m&T)] - 1 

in d = 3 (31) 



C(3/2) ,„ ^3/2 



This behavior determines the crossover phase boundaries 
shown in Fig^ The physical properties are determined 
by the larger of the two "mass" terms in the b Green's 
function, \fj,f,\ or £(,(0,0) - consequently, the crossover 



phase boundaries in Fig. 2] lie at T 



2/3 



Jk - 



Jkc\ 2 ' 3 - These boundaries separate the U(l) FL* region 
at low T and fib < 0, and the FL region at low T and 
fib > 0) from the intermediate quantum critical region. 
Note that there is no phase transition in the FL region at 
T > 0: this is due to the compactness of the underlying 
U(l) gauge theory, and the fact that the "Higgs" and 
"confining" phases are smoothly connected in a compact 
U(l) gauge theory in three total dimensions |2l| . 

We now briefly comment on the nature of the elec- 
trical transport in the three regions of Fig The be- 
havior is quite complicated, and we will first highlight 
the main results by simple estimates in the present sub- 
section. A more complete presentation based upon the 
quantum Boltzmann equation appears in Section IV CI 

The conventional FL region is the simplest, with the 
usual T 2 dependence of the resistivity — the gauge fluc- 
tuations are quenched by the "Meissner effect" . 

In the U(l) FL* region, there is an exponentially small 
density of thermally excited b quanta, and so the boson 
conductivity at is also exponentially small. As in earlier 
work |43j, the resistances of the b and / quanta add in 
series, and so the total b and / conductivity remains ex- 
ponentially small. The physical conductivity is therefore 
dominated by that of the c fermions, which again has a 
conventional T 2 dependence. 

Finally, we comment on the transport in the quan- 
tum critical region. This we will estimate following the 



method of Ref. with a more complete calculation ap- 
pearing in the following subsection. A standard Fermi's 
Golden rule computation of scattering off low-energy 
gauge fluctuations shows that a boson of energy e has 
a transport scattering rate 



Tbtr(e) 



(32) 



for energies e <C T 2 / 3 . From this, we may obtain the 
boson conductivity by inserting in the expression 



Ob 



d 3 fcr bt r(ebfc)fc 2 



dn(e bk ) 
de 



(33) 



where n(e) = l/(e e / T — 1) is the Bose function, and 
e bk = k 2 /{2m b ) + E 6 (0, 0) = k 2 /(2m b ) + c 2 T 3 / 2 for some 
constant c 2 . Estimating the integral in (|33|l we find that 
there is an incipient logarithmic divergence at small k 
which is cutoff by £(,(0, 0) ~ T 3 / 2 , and so 0& diverges 
logarithmically with T: 



a h ~ ln(l/T). 



(34) 



There are no changes to the estimate of the / conductiv- 
ity from earlier work H2,[2j|, and we have <7/ ~ T~ 5 / 3 . 
Using again the composition rule of Ref. I43I we see that 
the asymptoptic low-temperature physical conductivity 
is dominated by the behavior in l)34[l. 

As an aside, we note that for the theory H24|) in two spa- 
tial dimensions the result of Eq. I|34|l continues to hold, 
whereas the fermion part becomes 07 ~ J 1 - 4 / 3 . This im- 
plies that the asymptoptic low-T physical conductivity is 
dominated by (|34|l in d = 2 as well. 



C. Quantum Boltzmann equation 

We now address electrical transport properties of the 
theory (|24(l in more detail, using a quantum Boltzmann 
equation. The analysis is in the same spirit as the work of 
Ref. IH but, as we have discussed in Section lTAl the vari- 
ation in the boson density as a function of temperature 
leads to very different physical properties, and requires a 
distinct analysis of the transport equation. 

We saw in Section IV Bl that the electrical conductiv- 
ity was dominated by the b boson contribution, and so 
we focus on the time (t) dependence described by the 
distribution function 



(bl(t)b k (t)} 



(35) 



In the absence of an external (physical) electric field E, 
we have the steady state value f(k,t) = fo(k) with 



fo(k) 



(36) 



exp [(fc 2 /(2m b ) - fi b + £ h (0, 0))/T] - 1 ' 

with £(,(0,0) given in (|31|l . The transport equation in 
the presence of a non-zero E(t) can be derived by stan- 
dard means, and most simply by an application of Fermi's 
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FIG. 6: Plot of the function k 4 tp(k) for a few values of the 
reduced temperature T and the interaction parameter u 1401 . 

is defined in Eqs. 138H and I41L and has been obtained 
from the numerical solution of the quantum Boltzmann equa- 
tion 
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golden rule. The bosons are assumed to scatter off a fluc- 
tuating gauge field with a propagator given by Ijl7(l or 
()30J) . and this yields the equation 

df(k,t) p df(k,t) 



x(2ir)S 




f(k,t)(l + f(k + q,t))(l + n(n)) 



f(k + q,t)(l + f(k,t))n(Q) 



(37) 



where n(Q) is the Bose function at a temperature T as 
above. 

We will now present a complete numerical solution of 
(|37|l for the case of a weak, static electric field, to lin- 
ear order in E. The analysis near the quantum critical 
point parallels that of Ref. l45l with the main change be- 
ing that instead of the critical scattering appearing from 
the boson self-interaction u, the dominant scattering is 
from the gauge field fluctuations (note, however, that it 
is essential to include the interaction u to first order in 
the self-energy shift in Q36|l '). We write 



f(k,t)=f (k) + k-Ef 1 (k), 



(38) 



where notice that f\ depends only on the modulus of k 
and is independent of t. We now have to insert (|38|l into 
the transport equation (|37|l and the expression for the 
electrical current 



At) 



d d k k 
(2n) d m b 



mt), 



(39) 



FIG. 7: Scaling function for the boson conductivity, a = 
&b/(mbXf), as function of the reduced temperature T for dif- 
ferent values of the interaction parameter u 14011 . The re- 
sults are obtained from the numerical solution of the quantum 
Boltzmann equation 142H together with 1431 1. Top panel: con- 
ductivity a(T) on a log-log scale. Bottom panel: resistivity 
1/<t(T) on a linear scale. 



linearize everything in E, and so determine the propor- 
tionality between J and E. 

It is useful to re-write the equations in dimensionless 
quantities fl = f2/T, k = k/y/2m,bT, a = cri,/(mbXf)- 
Then it is easy to see that the solution of the quantum 
Boltzmann equation at the critical coupling, fif, — 0, is 
characterized by two parameters, 

f=^{2m b fT, 



COT r 

4tt 3 /2 



Xf 



(40) 



where T is a reduced temperature, and u parametrizes 
the temperature dependence of the effective "mass" of 
the bosons from Eq. I|31() : T and Xf are the parameters 
of the gauge propagator l|17|l . The linearized form of the 
Boltzmann equation (|37|l for the function 



h{k) = i>(k/^/2m b T) (41) 

is obtained as 

- f' a {k) = / dh [Kifc hMk) + K 2 (k, hMh)] (42) 
Jo 

with /o(x) = d / (dx 2 )[exp(x 2 +uVT) - 1] _1 ; the expres- 
sions for the functions K± t 2 are given in Appendix [El 
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From the solution of Eq. i|42|) one obtains the conductiv- 
ity according to 

1 f°° 

a(T,u) = =/ dkk A Mk). (43) 

The integral equation l|42ll was solved by straightfor- 
ward numerical iteration on a logarithmic momentum 
grid. We show sample solutions for the function k A tp{k) 
in Fig. The final results for the scaling function 
of the conductivity are displayed in Fig. [7| For small 
temperatures, the logarithmic divergence of o~b(T) an- 
nounced in Eq. l|34|l is clearly seen; for larger tempera- 
tures the conductivity is exponentially suppressed due to 
the temperature-dependent boson mass. In the crossover 
region, the results could be fitted with a power law over 
a restricted temperature range of roughly one decade, 
however, no extended power-law regime emerges. In 
comparison with experiments, one has to keep in mind 
that the physical resistivity is given by a sum of boson 
and fermion resistivities, and that the logarithmically de- 
creasing low-temperature part of l/<7b(T) cannot be eas- 
ily distinguished from a residual resistivity arising from 
impurities. 



D. SDW order 

Our discussion so far has focused primarily on the 
crossover between the FL to U(l) FL* phases, as this 
captures the primary physics of the Fermi volume chang- 
ing transition. At low T, we discussed in Section IIV Bl 
that the longitudinal part of the gauge fluctuations may 
induce SDW order on the spinon Fermi surface of the 
FL* phase (leading to the SDW* phase). On the FL 
side of the transition, the gauge fluctuations are formally 
gapped by the Anderson-Higgs mechanism. They will, 
however, still mediate a repulsive (though finite ranged) 
interaction between the quasiparticles at the hot Fermi 
surface. Furthermore the shape of the hot Fermi sur- 
face evolves smoothly from the spinon Fermi surface if 
the FL* phase. Consequently, it is to be expected that 
the SDW order will continue into the FL region up to 
some distance away from the transition. Thus it seems 
unlikely that there will be a direct transition from SDW* 
to FL at zero temperature. The actual situation then 
has some similarities to the mean-field phase diagram in 
FigEl However, fluctuations will strongly modify the po- 
sitions of the phase boundaries, and we expect that the 
U(l) FL* region actually occupies a larger portion of the 
phase diagram. Also there is no sharp transition between 
the FL and U(l) FL* regions (unlike the mean-field sit- 
uation in Fig0J, and there is instead expected to be a 
large intermediate quantum-critical region as shown in 
Fig. 13 



VI. EXPERIMENTAL PROBES OF THE U(l) 
SDW* STATE 

In this Section, we discuss experimental signatures of 
the U(l) SDW* phase focusing particularly on the dis- 
tinctions with more conventional SDW metals. 

We begin by considering a U(l) SDW* phase in which 
a portion of the spinon Fermi surface remains intact. As 
discussed in Section llVl the coupling between the gap- 
less spinons and the gauge field leads to singularities in 
the low-temperature thermodynamics in this phase. In 
particular the specific heat behaves as C(T) ~ Tln(l/T) 
at low temperature. Thus this phase is readily distin- 
guished experimentally from a conventional SDW. Elec- 
trical transport in this U(l) SDW* phase will be through 
the conduction electrons with no participation from the 
spinons. Thus electrical transport will be Fermi liquid- 
like. In contrast thermal transport will receive contribu- 
tions from both the conduction electrons and the gapless 
spinons. Consequently the thermal conductivity will be 
in excess of that expected on the basis of the Wiedemann- 
Franz law with the free electron Lorenz number. 

The distinction with conventional SDW phases is much 
more subtle for U(l) SDW* phases where the spinons 
have a full gap. In this case, there is a propagating gap- 
less linear dispersing photon which is sharp. The pres- 
ence of these gapless photon excitations potentially pro- 
vides a direct experimental signature of this phase. It is 
extremely important to realize that the emergent gauge 
structure of a fractionalized phase is completely robust 
to all local perturbations, and is not to be confused with 
any modes associated with broken symmetries. Thus de- 
spite its gaplessness the photon is not a Goldstone mode. 
In fact, the gaplessness of the photon is protected even 
if there are small terms in the microscopic Hamiltonian 
that break global spin rotation invariance. Being gap- 
less with a linear dispersion, the photons will contribute 
a T 3 specific heat at low T which will add to similar 
contributions from the magnons and the phonons of the 
crystal lattice. In addition, the conduction electrons will 
contribute a linear T term. The phonon contribution is 
presumably easily subtracted out by a comparison be- 
tween the heavy Fermi liquid and magnetic phases. To 
disentangle the magnon and photon contributions, it may 
be useful to exploit the robustness of the photons to 
perturbations. Thus for materials with an easy-plane 
anisotropy, application of an in-plane magnetic field will 
gap out the single magnon, but the photon will stay gap- 
less and will essentially be unaffected (at weak fields). 
Thus careful measurements of field-dependent specific 
heat may perhaps be useful in deciding whether the U(l) 
SDW* phase is realized. 

Finally, quasi-elastic Raman scattering has been sug- 
gested as a probe of the U(l) gauge field fluctuations 
|46j in the context of the cuprates — the same prediction 
applies essentially unchanged here to the fractionalized 
phases in d = 3. 

Conceptually the cleanest signature of the U(l) SDW* 
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phase would be detection of the gapped monopole. How- 
ever at present we do not know how this may be directly 
done in experiments. Designing such a "monopole detec- 
tion" experiment is an interesting open problem. 



VII. DISCUSSION 

The primary question which motivates this paper is 
how to reconcile a weak moment magnetic metal with 
non-Fermi liquid behavior close to the transition to the 
Fermi liquid. We have explored one concrete route to- 
ward such a reconciliation. The U(l) SDW* magnetic 
states discussed in this paper may be dubbed spin-charge 
separated spin density wave metals. They constitute a 
class distinct from both the conventional spin density 
wave metal and the local-moment metal mentioned in the 
Introduction. However, they share a number of similari- 
ties with both conventional metals. Just as in the conven- 
tional local-moment metal, in the U(l) SDW* state the 
local moments do not participate in the Fermi surface. 
Despite this the ordering moment may be very small. 
Indeed this state may be viewed as a spin density wave 
that has formed out of a parent non-magnetic metallic 
state with a "small Fermi surface" . This parent state is 
a fractionalizcd Fermi liquid in which the local moments 
have settled into a spin liquid and essentially decoupled 
from the conduction electrons. The spinons of the spin 
liquid form a Fermi surface which undergoes the SDW 
transition - this transition does not affect the deconfinc- 
ment property of the gauge field, because the SDW order 
parameter is gauge neutral and thus does not effectively 
couple to the gauge field excitations. 

We showed that in the region of evolution from this 
state to the conventional Fermi liquid, non-Fermi liq- 
uid behavior obtains (at least at intermediate temper- 
atures). We also argued that the underlying transition 
that leads to this non-Fermi liquid physics is the Fermi 
volume changing transition from FL to FL* . Despite the 
jump in the Fermi volume, this transition is continuous 
and characterized by the vanishing of the quasiparticle 
residue Z on an entire sheet of the Fermi surface (the 
"hot" Fermi surface) on approaching the transition from 
the FL side. 

A specific heat that behaves as Tln(l/T) is commonly 
observed in a variety of heavy-fermion materials close to 
the transition to magnetism. In the context of the ideas 
explored in this paper, such behavior of the specific heat 
is naturally obtained in three-dimensional systems. A 
small number of heavy-fermion materials exhibit such a 
singular specific heat even in the presence of long-ranged 
magnetic order. As we have emphasized, precisely such 
non-Fermi liquid specific heat obtains in one of the ex- 
otic magnetic metals discussed in this paper (the U(l) 
SDW* phase with a partially gapped spinon Fermi sur- 
face). It would be interesting to check for violations of 
the Wiedemann-Franz law at low temperature in such 
materials. 



A general point emphasized in this paper is that the 
observed non- Fermi liquid physics near the onset of mag- 
netism actually has little to do with fluctuations of the 
magnetic order parameter. Rather we propose that the 
non-Fermi liquid physics is associated with the destruc- 
tion of the large Fermi surface. The concrete realization 
of this picture explored in this paper is that the destruc- 
tion of the large Fermi surface leads to a fractionalized 
Fermi liquid which eventually (at low temperature) de- 
velops spin density wave order. As we discussed exten- 
sively, the resulting spin density wave state is an exotic 
magnetic metal. 

It is also of interest to consider a different scenario 
in which the small Fermi surface state is unstable at low 
temperature toward confinement of spinons and magnetic 
order. It is particularly interesting to consider such a sce- 
nario in d — 2. The physics of the Fermi- volume changing 
fluctuations is again described by a theory of condensa- 
tion of a slave boson field coupled to a Fermi surface of 
spinons by a U(l) gauge field. For a non-compact U(l) 
gauge field, such a theory has a number of interesting 
properties. As noticed by Altshuler et al. the spin 
susceptibility at extremal wavevectors of the spinon fermi 
surface have (possibly divergent) singularities due to the 
gauge fluctuations. Indeed the spin physics of this model 
is critical and described by a non-trivial fixed point. The 
dynamical susceptibility at these extremal wavevectors 
and at a frequency u> is expected to satisfy u>/T scal- 
ing. For a general spinon Fermi surface these extremal 
wavevectors will chart out one-dimensional lines in the 
Brillouin zone at which critical scattering will be seen in 
inelastic neutron scattering. A spin density wave insta- 
bility can develop out of this critical state at a particular 
extremal wavevector where the amplitude of the diverg- 
ing susceptibility is the largest. Arguments very similar 
to those in Section IV Bl also show that transport will be 
governed by non-Fermi liquid power laws in this theory. 

There is a striking qualitative resemblance between 
these results and the experiments on CeCue-^Auj;. At 
the critical Au concentration neutron scattering experi- 
ments see critical scattering on lines in the Brillouin zone 
satisfying uo/T scaling Magnetic ordering occurs at 
particular wavevectors on this line. Furthermore, em- 
pirically the spin fluctuations appear to be quasi two- 
dimensional, suggesting that the ideas sketched above 
may indeed be relevant. We note that they significantly 
differ from earlier proposals to explain the behavior of 
CeCu6-a;Au x |4^. As mentioned in the text, the spe- 
cific heat in the d = 2 quantum critical region will have 
the form C/T ~ T -1 / 3 ; interestingly such a behavior 
has been observed in YbRh2Si2 in the low-temperature 
regime near a quantum-critical point 49]. On the theo- 
retical front, there are a number of conceptual issues 
related to the legitimacy of ignoring the compactness of 
the gauge field in d = 2. Developing a more concrete the- 
oretical description of these general ideas is an interesting 
challenge for future work. 
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APPENDIX A: OSHIKAWA'S ARGUMENT AND 
TOPOLOGICAL ORDER 

Oshikawa has presented an elegant non- 

perturbative argument demonstrating that the volume 
of the Fermi surface is determined by the total number 
of electrons in the system. In our previous work p"H , and 
in the present paper, we have argued for the existence of 
a non-magnetic FL* state with a different Fermi surface 
volume. As we discussed earlier this apparent con- 
flict is resolved when we allow for global topological ex- 
citations in Oshikawa's analysis; such excitations emerge 
naturally in the gauge theories we have discussed for the 
FL* state. In other words, Oshikawa's argument implies 
that violation of Luttinger's theorem must be accompa- 
nied by topological order. 

In this Appendix, we briefly recall the steps in Os- 
hikawa's argument, and show how it can be modified to 
allow for a small Fermi surface in a FL* state. As far as 
possible, we follow the notation of Oshikawa's paper [5ll |. 

For definiteness, consider a two-dimensional Kondo 
lattice with a unit cell of lengths a x ,y The ground state 
is assumed to be non-magnetic, with equal numbers of up 
and down spin electrons. Place it on a torus of lengths 
L x ,y, with L x /a x , L y /a y co-prime integers. Adiabati- 
cally insert a flux <f> = 2ir (h = c = e = 1) into one of the 
holes of the torus (say the one enclosing the x circum- 
ference), acting only on the up-spin electrons. Then the 
initial and final Hamiltonians are related by a unitary 
transformation generated by 

^^(^E'm) ( ai ) 

where n r | is the number operator of all electrons (includ- 
ing the local moments) with spin up on the site r. After 
performing the unitary transformation to make the final 
Hamiltonian equivalent to the initial Hamiltonian, the 
final and initial states are found to have a total crystal 



momentum which differs by 

2tt L x L y p a ( ,271- \ 

AP X = y -— mod— (A2) 

L x v 2 \ a x ) 

where vq — a x a y is a volume of a unit cell, the second 
factor on the r.h.s. counts the number of unit cells in the 
system, and p a = 2p a j is the mean number of electrons in 
every unit cell. Clearly the crystal momentum is defined 
modulo 2ir/a x , and hence the modulus in (|A2|) . 

Now imagine computing the change in crystal momen- 
tum by studying the response of the quasiparticles to the 
inserted flux. As shown by Oshikawa, the quasiparticles 
associated with a Fermi surface of volume V lead to a 
change in momentum which is 

Ap -' - 1 wTO ( mo ^ ) • (A3) 

where the second factor on the r.h.s. counts the number 
of quasiparticles within the Fermi surface. Equating AP X 
and AP£, and the corresponding expressions for AP y 
and AP^, Oshikawa obtained the conventional Luttinger 
theorem, which applies to the volume V = Vfl of the 
Fermi surface in the FL state 

2 (2^ VFL = Pa(m ° d 2) (M) 

In the FL* state, there are additional low-energy ex- 
citations of the local moments that yield an additional 
topological contribution to the change in crystal momen- 
tum. Indeed, the influence of an insertion of flux $ 
is closely analogous to the transformation in the Lieb- 
Schultz-Mattis [52] argument, and it was shown |53l 
that a spin liquid state in d = 2 acquires the momentum 
change 

AP* = 1-h. Lod— ^) (A5) 
a x a y \ a x J 

where the second factor on the r.h.s. now counts the 
number of rows which have undergone the Lieb-Schultz- 
Mattis transformation. Now using AP X = APf + APj, 
we now obtain the modified Luttinger theorem obeyed in 
the FL* phase: 

2 7T^2 = - l)(mod 2). (A6) 

It is clear that the above argument is easily extended 
to a Z 2 FL* state in d = 3. The case of U(l) FL* state 
in d = 3 is somewhat more delicate because there is now 
a gapless spectrum of gauge fluctuations which can con- 
tribute to the evolution of the wavefunction under the 
flux insertion; nevertheless, the momentum change in 
(|A5|) corresponds to an allowed gauge flux, and we ex- 
pect that (|A5|) continues to apply. 
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APPENDIX B: TOY MODEL WITH U(l) 
FRACTIONALIZATION AND A SPINON FERMI 
SURFACE 

In this Appendix, we will display a concrete model in 
three spatial dimensions that is in a U(l) fractionalized 
phase in three dimensions, and has a Fermi surface of 
spinons coupled to a gapless U(l) gauge field. As dis- 
cussed earlier, this spinon Fermi surface could eventually 
(at low energies) undergo various instabilities including 
in particular to a spin density wave state. 

Consider the following model: 

H = H ti ,+H A +H b + H u + Hu, (Bl) 

H ti f, = - t (y4vv + h.c.) , 

(rr') 

H A = A e ^"' (f^l - 4^ti) + h-c. , 

(rr') 

Hb = — W COs((f) rr f — (j>rr ,f ) 7 

[rr'r"] 

H u = u } j n 2 r , , 

(rr 1 ) 
r 

Here ipr destroys a spinful charge- 1 electron at each site of 
a cubic lattice in three spatial dimensions; e lc ^ rr ' creates a 
charge-2, spin-0 "Cooper pair" that resides on the links. 
n rr i is conjugate to c/) rr i and may be regarded as the 
Cooper pair number associated with each link. N r is the 
total charge associated with each site and is given by 

N r = ^2 nrr ' + ^Vv- (B2) 

The Hamiltonian H may be regarded as describing a sys- 
tem of electrons coupled with strong phase fluctuations. 
The first term in Hb represents Josephson coupling be- 
tween two "nearest neighbor" bonds. H u penalizes fluc- 
tuations in the Cooper pair number at each bond. Hjj 
penalizes fluctuations in the total charge N r that can be 
associated with each lattice site. The total charge of the 
full system clearly is 

N tot =J2 N r- (B3) 

r 

Depending on the various model parameters, several 
distinct phases are possible. Here we focus on the limit 
of large U. Diagonalizing Hjj requires that the ground 
state(s) satisfy N r = 1 at all sites r. There is a gap 
of order U to states that do not satisfy this condition. 
Clearly the system is insulating in this limit. 

The condition N r = 1 for all r still allows for a huge 
degeneracy of ground states which will be split once the 
other terms in the Hamiltonian are included. This split- 
ting may be described by deriving an effective Hamilto- 
nian that lives in the space of degenerate states specified 



by N r — 1 . As discussed in Ref . IT3 . this effective Hamil- 
tonian may be usefully viewed as a (compact) U(l) gauge 
theory. This may be explicitly brought out in the present 
case by the change of variables 

n rr i = e r E rr i (B4) 

Ipra = fra for T S A , 

A* = Kfsflp fOT r e B . (B5) 

Here e r = +1 on the A sublattice and —1 on the B 
sublattice. In terms of these variables, the constraint 
N r = 1 reads 

V-E + flf r = l (B6) 

at each site r. (We note that both a and E may be 
regarded as vector fields defined on the lattice). At order 
w 2 /U, u, t, A, the effective Hamiltonian takes the form 

-Heff = Hk + H u + Hf , (B7) 

H K = -if J^cos(V x a) , 
p 

r 

H f = -A^2 (e ia -'/r + /r<+h.c.) . 

(rr') 

Here K = 2w 2 /U, and the sum ^ p runs over elemen- 
tary plaquettes. H c g together with the constraint may be 
viewed as a Hamiltonian for a compact U(l) gauge the- 
ory coupled to a gauge charge-1 fermionic matter field 
/. (Note that to leading order the t term does not con- 
tribute). H c ff still admits several different phases de- 
pending on its parameters. Of interest to us is the limit 
K ~ w 2 /U ^> u. In this limit, monopoles of the compact 
U(l) gauge field will be gapped. Consequently at low 
energies, we may take the gauge field to be non-compact. 
The cos(V x a) term can then be expanded to quadratic 
order to get the usual Maxwell dynamics for the gauge 
field. The /-particles form a Fermi surface which is cou- 
pled to this gapless U(l) gauge field. Note that in the 
low-energy manifold with N r — 1 at all r, all excitations 
have zero physical electric charge. Thus the / particles 
are neutral fermionic spinons. 

As with any Fermi surface, this spinon Fermi surface 
state could at low energies further undergo various insta- 
bilities to other states (density waves, pairing, etc) de- 
pending on the residual interactions between the spinons. 
There are various sources of such interactions: First there 
is the gauge interaction that is explicit in H c g in the lead- 
ing order. Then the t term contributes to H c s at second 
order and leads to a quartic spinon-spinon interaction as 
well. The specific low-energy instability of the spinon 
Fermi surface will be determined by the details of the 
competition between these various sources of interaction, 
and will not be discussed further here for this model. 
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FIG. 8: Mean-field phase diagram of H m t (|18|l on the cubic 
lattice, now as function of Kondo coupling Jk and external 
field H z at T = 0. Parameter values are as in Fig. 2] For a 
description of the phases see text. 



Apart from the deconfined phase discussed above, the 
model possesses confined phases; for large U those occur 
for smaller u, and the deconfinement transition occurs 
through the condensation of monopoles in the gauge field. 



APPENDIX C: MEAN-FIELD PHASE DIAGRAM 
IN AN EXTERNAL ZEEMAN MAGNETIC FIELD 

In this Appendix we briefly discuss the behavior of the 
mean-field theory of Sec. IIV CI in an externally applied 
field. A sample zero-temperature phase diagram is dis- 
played in Fig. |SJ which shows very rich behavior. 

The phases at small fields are straightforward general- 
izations of the low-temperature zero-field phases of Fig. 0] 
The U(l) SDW* has weakly polarized conduction elec- 
trons, bo = 0, non-zero xo indicating spinon hopping, and 
a canted spinon magnetization M r with a staggered com- 
ponent along x and a uniform component along z. The 
SDW phase has similar characteristics, but now bo ^ 
indicating a conventional weakly field-polarized magnet 
with confinement. Finally, the FL phase has bo ^ 0, 
Xo 7^ 0, weakly polarized heavy quasiparticles, and the 
mean-field parameter M r has only a uniform z compo- 
nent. 

In the small- Jk region, increasing external field pro- 
gressively suppresses the effect of Jh- At intermediate 
field, a phase with "canted" / moments arises, where now 
Xo = bo = (no spinon hopping), and M r is canted as 
described above. Larger fields fully polarize the local mo- 
ments, i.e., M r points uniformly along z with maximum 
amplitude, and Xo = &o = 0. This phase is also realized 
for larger Jk and large fields - here the field quenches the 
Kondo effect. On the Fermi-liquid side of the phase dia- 
gram, two more phases arise in the present mean-field 



theory which are labelled by FL2 and FL3 in Fig. [SJ 
both have non-zero bo and xo- In FL2, the mean- field 
parameter M r has both staggered and uniform compo- 
nents, i.e., this phase describes canted, weakly screened 
local moments. Turning to the FL3 phase, this high-field 
phase has the same symmetry characteristics as FL at 
intermediate fields, but a different Fermi surface topol- 
ogy. Whereas FL phase at intermediate fields has a single 
Fermi surface sheet for one spin direction (the major- 
ity spin have one full and one empty band whereas the 
minority spins have one partially filled and one empty 
band), in FL3 the upper band of the majority spins be- 
comes partially filled, too. FL and FL3 are separated 
by a strongly first-order transition in mean-field theory. 
There are numerous other phase transitions associated 
with a change in the Fermi surface topology - those do 
not display strong thermodynamic signatures and are not 
shown. We note that for the field range displayed in 
Fig- El \H cx t\ <C t, the conduction electrons are in gen- 
eral weakly affected by the field; significant polarization 
of them occurs only at much higher fields. 

Notably, smaller values of the decoupling parameter x 
admit yet another field-induced transition in the small- 
Jk region: If the magnetism is very weak, i.e., the 
spinons have a small gap compared to their bandwidth, 
then a small applied field can close the spinon gap with- 
out significantly affecting their band structure. Such a 
transition would yield a kink in the magnetization of the 
local-moment subsystem as function of the applied field, 
implying a "metamagnetic" behavior which is here gener- 
ically associated with a continuous transition. 



APPENDIX D: SPECIFIC HEAT SINGULARITY 

Here we present some details on the calculation of the 
singular specific heat coming from gauge fluctuations. 
The calculations in the FL* phase and at the critical 
point are standard. We will therefore only consider the 
FL phase. In this phase close to the critical point, trans- 
verse gauge fluctuations are described by the action 

*=/|^e(^+^><^..>i 2 - (di) 

As explained in Section^ close to the transition p s ~ 6g. 
This gives a free energy 

-|e/0.»(¥^+4 (D2 » 

To calculate the low-temperature specific heat, we need 
the change in free energy on going from zero to a small 
non-zero temperature. After a Poisson resummation this 
is given by 

8F(T) = F(T)-F{0) (D3) 
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= 2 



iftmuj— u( \u | +p s +A)) 



The uj, A integrals may now be performed to obtain 



equation 142|) in Sec. IV CI Inserting the ansatz (|38(l into 
Eq. Ij37|l leads to a scalar equation for fi . The frequency 
integral is easily performed; the remaining momentum 
integral can be split into radial and angular part. This 
directly yields Eq. , with 



K^kM) = [l + f {h)+n(k 2 -k 2 )] 



(El) 



a r j^l r P -uk(k +p s ) 

SF(T) = - Y du „ . — 



dkk 2 f (iruT coth{TruT) - \) e -^Hk 2 + Ps ) 



2tt 3 



In the last equation we have introduced an upper cut-off 
A for the momentum integral. The remaining integrals 
can now be straightforwardly evaluated for small T, and 
we find 



T 2 /A 2 

SF(T) = In — 

V ' 12tt \p s 



Thus the specific heat 



C(T) = 7 T 



(D4) 



(D5) 



with 7 ~ ln(l/p s ) ~ ln(l/6o)- Setting p s = 0, a similar 
calculation also shows that C(T) ~ Tln(l/T) in the FL* 
phase. 

For completeness, we mention the corresponding be- 
havior in two dimensions. In analogy to the above cal- 
culations, we find C{T) oc T 2 / 3 in the quantum-critical 
and FL* regions. 



K 2 (k,h 



[f {k)-n{k 2 -kl)} 



(E2) 



k k? f 1 , xki ,, / fci r / - 
with /o(x) = [exp(x 2 + uvT) — 1] 1 and n(x) = (e x 



1) . The kernel K(x, ki/k, kvT) is given by 



K(x, a, A) = Im —i 



l~a 2 



x 1- 



^/a 2 + T—2ax 
(ax- l) 2 



+ A(a 2 + l-2ax) 



a 2 + 1 — 2aa; 



The integrals necessary for the evaluation of Ki t 2 are of 
the form 



dy 



1-1/2 



y 3 + c 



with n = 0, . . . , 3 and can be performed analytically. 

The numerical solution of Eq. (|42|l is done by re- writing 
it in the form 



APPENDIX E: DETAILS OF THE QUANTUM 
BOLTZMANN EQUATION 

In the following we describe a few details of the deriva- 
tion of the linearized version of the quantum Boltzmann 



to 



+ K 2 (k,h) 



1>W 



which allows for a stable numerical iteration. 
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